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Abstract 

We investigate axially symmetric asymptotically flat vacuum self-gravitating system. 
A class of initial data with apparent horizon was numerically constructed. The examined 
solutions satisfy the Penrose inequality. The prior analysis of a massive system and the 
present results suggest that either massive or sourcefree configurations fulfill the Penrose 
inequality. 
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1 .Introduction. 

The problem of the Cosmic Censorship Hypothesis is one of the most significant issues 
of the classical gravity. The feasible technique of testing it is the examination of the validity 
of the Penrose inequality [1] . Breaking of the Penrose inequality points to the violation of 
the Cosmic Censorship. 

The geometry created by massive nonspherical (spheroidal) bodies has been numeri- 
cally analyzed in [2] . The Penrose inequality in this case has been satisfied. We would like 
to stress out that those calculations was made under the assumption of conformal flatness. 
Conformal flatness may be intuitively understood as a restriction to the case without grav- 
itational waves - a conformally flat system without matter is flat. Therefore the subject 
of our interest in this paper is in some sense an orthogonal situation: a geometry of the 
collapsing pure gravitational wave. After a few numerical tests [3] it has been recently 
analytically established [4] that vacuum solutions with apparent horizon exist. 

In this paper we are numerically examining the following formulation of the inequality: 

Let E be an asymptotically Eat Cauchy surface with a non-negative energy density, m 
- the Arnowitt-Deser-Misner mass and S the area of an external apparent horizon. Then 



Penrose argues that if the above is not true then the Cosmic Censorship will be broken. 
A more refined version of (1) has been proposed in [5]. The inequality (1) is still proven 
only for a restricted class of geometries [6,7] (see also a short review in [8]). 

2. Main equations. 

We examine the formation of apparent horizons on the initial hypersurface E which 
is required to be momentarily static and axially symmetric (Brill waves). In the case 
of vacuum constraints on E reduces to 




(1) 




Due to vanishing of the scalar curvature, $ and q are constrained by [9] 



A$ + = 0, 



(2) 



where A denotes 3-dimensional fiat laplacian and 




(3) 
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Regularity of metric on the axis implies the boundary conditions for $ and q. 



^{e = o) = o, |^(^ = o) = o, ^(^ = o) = o. (4) 

The equivalent group of conditions must hold for 9 = 71. Moreover asymptotically solution 
is flat: $^1 + ^,5-^0. Where m is the ADM mass (see equation (1)). The above 
requirements may be satisfied by following choices of function q [3]: 

qi{r,e)^ ^^^^ . (5) 

Inserting q in (3) we obtain 

/i = - ^ (2C/ - 45C/ sin^ 9 + 50C/ sin^ 9) , (6) 
8 

where 

Similarly we analyze another q: 

q2{r,9) = Ae-'^''\'^sm^9. (8) 



Now / has the following form: 

/2 = - je"'^^' (1 - Gar^ sin^ 9 + 2a sin^ 9). (9) 

The above formulas we use in the next analysis for calculation of function $. Because of 
the regularity the determinant of the metric must have no zeros. Therefore we consider 
only solutions which are everywhere positive; $ > 0. 

A surface S is called an apparent horizon if the expansion of outgoing future 
directed null geodesies which are orthogonal to S vanishes everywhere on S. In the mo- 
mentarily static case expansion is the divergence of normal unit vector to S. 

Din' = 0. (10) 



Hence on time symmetric S notion of apparent horizon coincides with definition of minimal 
surface. We investigate the external apparent horizon which is the outermost apparent 
horizon that surrounds the region of concentration of gravitational energy. Equation (10) 
for our metric takes the form 

r|/4$^ q^ \ 2/4$^ 
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+^^(^-y + coMj-r (^-y + -)=0 (11) 



The area of S is given by 



TT 

S^27r j $^e-«/2^r^ rsin« (12) 



3. Numerical techniques and results of calculations. 

Originally we intended to solve the eq. (2) by the use one of standard method on 
lattice. The numerical error was too big and therefore we solve the partial differential 
equation (2) by expanding the functions $ and / in a series of the Legcndrc polynomials 
[3]. We obtained an infinite system of ordinary differential equations for r dependent 
functions. We approximately solved this system by retaining only leading terms (up to 
1=5). Numerical results shows that the contribution from the higher terms is very small 
and safely can be neglected. Our numerical analysis of the system just described and the 
equation for external apparent horizon (11) was based of the fourth order Rungc-Kutta 
method. For q2 we examined only one value of a = 0.1. The solutions for other values of 
a may be obtained from the previous ones by simple scaling transformations. 

We analyzed only such configurations for which the shape of the external apparent 
horizon significantly diverges from spherical symmetry. In the case of spherical symmetry 
the inequality (1) reduces to an identity. In the below tables we collect the most important 
results of our calculations. The first table corresponds to q given by formula (5) and the 
second one to (8) with a = 0.1. 



A 


m 


ro 




S 


m — s/S/lQiT 
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4.59719 


2.13920 


2.08129 


1055.72607 


0.014 


8.5 


7.30543 


3.63536 


3.59755 


2681.55505 


0.0015 
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14.00554 


7.00300 


6.98652 


9858.96851 


0.0006 



table 1. 



A 


m 


ro 


rTT/2 


S 


m, — \J SjX^Ti 


-2.1 


11.76569 


5.13853 


7.90717 


10905.86182 


0.0359 


-2.5 


16.63314 


6.47847 


8.85901 


13864.95839 


0.0249 


-2.6 


18.77306 


7.85006 


9.86873 


17680.94401 


0.0180 



table 2 

Summarizing, the obtained result confirm the validity of the Penrose inequality for the 
analyzed class of Brill waves. Although the proof of the inequality (1) is still uncompleted 
it seems that our belief in this conjecture is rather reasonable. 
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